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We study pairwise quantum entanglement in Systems of fermions itinérant in a lattice from a 
second-quantized perspective. Entanglement in the grand-canonical ensemble is studied, both for 
energy eigenstates and for the thermal state. Relations between entanglement and superconducting 
corrélations are discussed in a BCS-like model and for 77-pair superconductivity. 
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PACS numbers: 03.67.Lx, 03.65.Fd 
I. INTRODUCTION 

The concept of quantum entanglement ^ is believed 
to play an essential rôle in quantum information pro- 
cessing (QIP) As a conséquence many efforts have 
been devoted to the characterization of entanglement . 
The very définition of entanglement relies on the ten- 
sor product structure of the state-space of a composite 
quantum System. Unfortunately, due to quantum statis- 
tics, such a structure does not appear in an obvions fash- 
ion for Systems of indistinguishable particles, i.e., bosons 
or fermions. Indeed for thèse Systems, in view of the 
(anti)symmetrization postulate, one has to restrict to 
a subspace of the A^-fold tensor product of the single 
particle spaces. Such a subspace, e.g., the totally anti- 
symmetric one, has not a naturally selected tensor prod- 
uct structure. It turns ont that the notion of entangle- 
ment is affected for Systems of indistinguishable particles 
by some ambiguity. 

Since it is of direct relevance to several implementa- 
tion proposais for QIP e.g., quantum-dots based, this 
issue has been very recently addressed in the literature 
[ ï . A quantum computation model was proposed 

[ 3| by using L local fermionic modes (LFMs)- sites which 
can be either empty or occupied by a fermion. Moreover 
the use of quantum statistics for some QIP protocols have 
been analyzed 0. 

Along the same line of realizing a bridge between quan- 
tum information science and conventional many-body 
theory it has been discussed entanglement in magnetic 
Systems [ pT| , p^ljr^JÏ^ , [Ï5|| . In particular entanglement in 
both the ground state ]ll| , p^ and thermal state ||ï^,^, Ï5| 
of a spin-1/2 Heisenberg spin chain have been analyzed in 
the literature. In this situation the System state is given 
by the Gibb's density operator pT = exp{—H/kT) /Z, 
where Z =tr[exp (— iï/feT)] is the partition function, H 
the System Hamiltonian, k is Boltzmann's constant which 
we henceforth will take equal to 1, and T the tempéra- 
ture. As p(T) represents a thermal state, the entangle- 
ment in the state is called thermal entanglement [ [13| . 
Finally the intriguing issue of the relation between en- 
tanglement and quantum phase transition |16| have been 
addressed in a few quite récent papers p^ , p8| . 

In this paper we will explore the relations between en- 
tanglement and (super)conducting corrélations by foUow- 
ing the spirit of Réf. ||]. It is important to stress that. 



due to the lack of measure of genuine many-body en- 
tanglement, we restrict us pairwise entanglement in this 
paper. Notice that in this approach the subsystems are 
given by modes and not by particles. This is therefore an 
essentially second-quantized appTO&ch [|T]. 

In Sec. II basic définitions are given and the mapping 
scheme between LFMs and qubits introduced in ^ is 
briefiy recalleed. In Sec. III the entanglement in both 
eigenstates and thermal state is studied for free fermions 
hopping in a lattice. In Sec. IV the relations between 
pairwise entanglement of and superconducting corréla- 
tions are discussed for Two typ es of superconductivity, 
BCS-like superconductivity |l9|, and the so-called 77-pair 
superconductivity [pÔl. Sec. V contains the conclusions. 



II. LATTICE FERMIONS 

Let us start by recalling some basic facts about (spin- 
less) fermions on a lattice. In the second-quantized pic- 
ture the basic objects are the création and annihilation 
operators cj and Ci of l-th LFM. They satisfy the canon- 
ical anti-commutation relations 



[ci,Cj]+ = 0, [ci,c]]+ = %. 



(1) 



The Hilbert space naturally associated to the L LFMs, 
known as Fock space Hf, spanned by 2^ basis vectors 
|ni,...,ni) :=nti(4)"'|0) {ni=0,l\fl). 

From the above occupation-number basis it should be 
évident that Tip is isomorphic to the L-qubit space. This 
is easily seen by defining the mapping 



A:=f[i4r\0) 



^1 



ni) 



i(^;+no), (2) 



1=1 



where is the raising operator of l-th qubit. This 
is an Hilbert-space isomorphism between 7ip and C^^. 
By means of this identification one can discuss entangle- 
ment of fermions by studying the entanglement of qubits. 
Clearly this entanglement is strongly relative to the map- 
ping (g) and it is by no means unique. By defining new 
fermionic modes by automorphisms of the algebra defined 
by Eq.(^ one gives rise to différent mappings between 
Hf and C®'^ with an associated différent entanglement. 
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This simple fact is one of the manifestations of the rela- 
tivity of entanglement . 

It is useful to see how the mapping (^) looks on the 
operator algebra level. From the relation c||ni, til) = 



^r^i - 1 

'5«i,o(-l)^'==i' 
that 



= |ni, ...,ri;_i, 1, 7i;+i, ...,nL), it foUows 



i-i 



(3) 



fe=i 



where is the z component of the usual Pauli matri- 
ces for fc-th qubit. This algebra isomorphism is quitc 
well-known in the condensed matter literature and it 
referred to as the Jordan- Wigner transformation [ p3[ . 
Notice that the inverse of Eq.(H) is given by aj t-^ 

c[ J||j,^\exp(i7rc^Cfc). We see that, due to the non-local 
character of the mapping A(A^^) even simple fermionic 
(spin) models can be transformed into non-trivial spin 
(fermionic) models. On the other hand the fermionic 
state likc Y[k ^IlO) are clearly mapped by A onto product 
qubit States. 

It is important to keep in mind that, for charged 
and/or massive fermions, the Fock space is not the state- 
space of any physical System. Indeed, at variance with 
massless neutral particles, e.g., photons, only eigenstates 
of iV = X^i^i are allowed physical vectors and, for the 
same reason, only operators commuting with N could be 
physical observables. This of course is nothing but a su- 
perselection rule, i.e., TCp = QN^oT^piN) that does not 
allow for linear superposition of states corresponding to 
différent charge/mass eigenvalues |Q. 

Despite the above considérations we notice that in 
some situations one is led to attribute to the whole Fock 
space some physical meaning. This happens for Systems 
in a symmetry broken phase. For example in super- 
conductivity and superfuidity the order parameter corre- 
sponds to an expectation value of an operator Connecting 
différent iV-sectors. It follows that the associated mean- 
field Hamiltonian does not commute with N. 

Of course one can argue that this kind of violation oc- 
curs on a level that has not any deep physical significance, 
after ail the mean-field approach is just a variational one 
aimed to produce good approximation to physical expec- 
tation values. According to this view therefore the prop- 
erties, e.g., entanglement, of the ansatz states should not 
to be taken too seriously. Nevertheless we think that this 
issue has some interest and the relations between pairwise 
entanglement and superconductivity will be provided be- 
fore the section of conclusions. 



III. ITINERANT SYSTEMS 

Let us now consider free spinless fermions in a lattice. 
The Hamiltonian is given by 



H ^-tJ2 (c^Q+i + 4+iCi) -l^J2 



1=1 



(4) 



1=1 



with the periodic boundary condition. Here t represents 
the hopping intégral between sites and /i is the chemical 
potential. 

It is known that the eigenvalue problem of H can be 
solved by a discrète Fourier transformation (DFT) 



(5) 



fc=i 



where ui — exp(i27r/L). After the DFT, the Hamiltonian 
(0) becomes 



H = -2tY^ cos(27rfc/L)4cfc - fiN, 



(6) 



k=l 



where N — ^I^fc ~ cjci is the total ferniion 

number operator. From Eq.(|6|), we immediately obtain 
the eigenvectors 



kjv = {ki,k2, ...,kN) e Z^, 
and the corresponding eigenvalues 



N 



EkN = ^i'^k, - Ai), Éfc, = -2tcos{2nki/L). 



(7) 



(8) 



1=1 



Associated to the new fermionic modes Ck there is a 
tensor product structure for the Fock space. The latter 
is defined by the mapping 



Adft :-n(4n0> 



(9) 



k=l 



Obviously since the eigenstates |kAr) are products with 
respect the tensor product structure due to Adft the en- 
tanglement in the eigenstates (0) is always zéro. However 
entanglement, associated with map A, may exist in the 
eigenstates. For instance, the concurrence C = 2/L for 
any pair of fermions when N = 1. The corresponding 
eigenstates are called W states pa,M,|4[ . 



A. Entanglement in the eigenstates 

In order to make an analysis of the entanglement of our 
spinless fermions we will use the notion of concurrence 
p7| , This is a simple measure for two qubits that allows to 
quantify the entanglement between any pair of fermions 
by our mapping. 

We define the reduced density matrix associated to the 
first and second LFMs as p^^^' eEnd(C^). Note that 



2 



the Hamiltonian is translation invariant, therefore entan- 
glenient between nearest-neighbor fermions are identical. 
Due to the fact that [N,H] — 0, the redueed density 
matrix have the foUowing form 



P 



(12) _ 



/ 



Wi z 

Z* W2 



(10) 



The nonvanishing relevant matrix éléments of p^^"^^ are 
given by ((•) dénotes the expectation value over p), 

u^l-2{N)/L+{nin2), 
V = {nin2), 

Z = {C\C2). 



The concurrence of p'^' is then given by 
C = 2 max{0, \z\ — ^/mj}. 



(11) 



(12) 



As the matrix éléments wi and W2 do not appear in the 
concurrence (p^), their expressions are not given in this 
paper. From Eqs.(|ïï|) and ( |Ï2|) it foUows that, in or- 
der to obtain the concurrence, we need to compute the 
corrélation functions {nin2) and {c\c2) and the mean 
fermionic number (TV). For the eigenstate jk^r) , after 
direct calculations |2al, we obtain 



u 

V 

z 



(n- 1)2-15 



tJV I I 



where n — N/L is the filling and 



N 

1=1 



N 



^iki 27r / L 



(13) 



(14) 



1=1 



By combining Eq. (|12|) and (|13|) one gets the concur- 
rence between two LFMs, 



C = 2max{0,|5k„|- 

{[{n~l)'-\SuJ'][n'-\S^J^]y/'}, 



(15) 



which is determined only the filling and the corrélation 
function {c[c2) — S'kjv- This latter quantity is obviously 
related to the "itinerancy" of the state, i.e., how fermions 
propagate. It foUows that the concurrence ( p^ ) contains 
direct information about the conducting properties of the 
given quantum state. 

From Eq.(p^ one can directly see that there exists 
entanglement if the corrélation function z and the filling 



factor n satisfy the équation 1^1 — 2{n 



l)\z[ 



[ri^ — nY < 0. Then we obtain that there exists pairwise 
entanglement between two LFMs in the eigenstates ii\z\'^ 
in the range — n + 1 — \/2v? — 2n + 1 < |zp < r? — 
71 + 1 + ^J2v? - 2ri + 1. In the case of TV = 1, |S'k„| = 



n = 1/L, and therefore C = 2/L ||J|^,|Ï|. For L = 2, 
this state is maximally entangled. A trivial case is = 0, 
there is no entanglement at ail. Now we consider the case 
N = L, i.e, the lattice is fuUy filled. Now |S'kjvl = 0, and 
hence C = 0. For the half filling (n = 1/2), u — v and 
Eq.(|ï|) reduces to C = 2max{0, |5k„| + l^k^P - 1/4}. 
Then the entanglement exists if and only if {^/2 — l)/2 < 
|5k«|<l/2. 

Now we consider the ground-state \G) of the System. 
It is obtained by filling the lowest single particle en- 
ergy levels. By taking for simplicity A^ odd one has 



|G) 



rN/2- 



Un=l 



1 t 



|0)(fc„ 



2Tm/L). From this 

définition it follows Sq = L-^{1 + 23?^,^^?"^^") = 
2L-i[cos(7r(Af/2 - 1)/L) sin((Af/2 - 1))/ sin(7r/L) - 1], 
now we take the limit L >—>■ oo, N/L i—>- n, the resulting 
expression is given by 

Soin) — — sm{TTn) (16) 

TT 

where we have used limi^oo ^ sin(7r/L) = tt. Moreover 
from the inequality Sq > n (1 — n) it follows that the 
second argument of the max function in Eq. (??) it is 
always non négative and we can get rid of the maximiza- 
tion. the concurrence then becomes 

C = 2{sin(7rn)/7r- [(n - 1)2 -sin2(7rn)/7r2]i/2} 

x[n2-sin2(7rn)/7r2]i/2}^ (17) 

in the infinité lattice. Fig. 1 shows the concurrence as 
a function of the filling n in the infinité lattice. We see 
that the entanglement becomes maximal at half filling 
for two ncighboring fermions, and the entanglement is 
symmetric with respect to the point of half filling. At the 
point of half filling the concurrence simply becomes C = 
2/tt + 2/tt'^ - 1/2 « 0.339262. The property of symmetry 
can also be seen directly from Eq.(^^ as the concurrence 
is invariant if we make the transformation n i— > 1 — n. 

For zéro chemical potential i.e., half- filling , it exists a 
direct relation between the concurrence and the ground 
state energy density êq. Indeed from translational sym- 
metry of the hamiltonian one has {c\c2) = -l/2tL{H) = 
— l/2teo(n), where we have used even the reality of {c\c2- 
This relation, that can be extend to finite température as 
well, is in a sensé remarkable in that it connects entangle- 
ment with a thermodynamical quantity that dépends on 
just the partition function of System. The latter is deter- 
mined just by the Hamiltonian spectrum, whereas Com- 
puting concurrence requires in gênerai also the knowl- 
edge of the eigenstate. This kind of connection between 
entanglement and thermodynamical quantifies has been 
discussed even for spin chains, both zéro [?] and finite 
température [pÔl. 



B. Thermal Entanglement 

In this section we extend our analysis to the entangle- 
ment in a finite température. The state of fermions at 
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thermal equilibrium is described by the foUowing Gibbs 
gran-canonical state 



PT = ^^exp(-/3£;k«) |kjv)(kAr| 



(18) 



where (3 — l/kT, k is the Boltzman's constant, and the 
partition function Z is given by 

L 

Z = ^exp(-/3£;k„) = Hil + e-f'^^''-^^ (19) 

kiv k—1 

The average accupation numbers (rife) is given by 



{nk)T 



1 



^Pick-p.) _|_ l ' 



(20) 



which is the Fermi-Dirac distribution. The expecation 
value (iV) is then easily obtained as {N) = J2k=ii^k)- 

(12) 

The reduced density matrix pj, associated with Adft 
will be a 4 X 4 diagonal matrix. Then from Eq.(p^), the 
concurrence is zéro for any pair of LFMs. So we will 
discuss the entanglement associated to the map A. The 
form of the reduced density matrix in the state pT is then 
given by Eq.(|ÏÔ|) and the concurrence is given by Eq.(12). 
Now we need to calculate the corrélation functions in Eq. 
( |rï| ) for the state pr- The resuit is readily obtained from 
Eq. ( p^ ) by replacing the eigenvalues of the n^'s with 
the corresponding thermal averages: 



{c\c,)=L-'J2'''(''k). 



k=l 



{n,n2) = {Nf/L'-\{c\c2)\' 



(21) 



Then the exact expression for the concurrence is given 
by the combination of Eqs.([ïï]), (|Ï2|), and Note the 
concurrence is thus obtained in an analytical form for 
arbitrary N. To exemplify this resuit let us first consider 
the simple case of L = 2 For two sites, from Eqs.(pT|), 
(p^, and (PU), the concurrence is given by 



C 



max{0,sinh(2/3|t|) - 1} 
cosh{f3p) + cosh(2/3i) ' 



(22) 



which is similar to the concurrence in a thermal state of 
the two-qubit XX model [Q. For three spécial values 
of (N), we have table I, from which we see that there is 
no entanglement when the chemical potential is =Foo (the 
corresponding mean fermion number is and 2). The en- 
tanglement is maximal when p = and other parameters 
are fixed. The mean fermion number is simply given by 



e''^ + cosh(2/3i) 
cosh{l3p) + cosh(2/3i) ' 



(23) 



from 

which we obtain ^ = /J-i ln{(2-(7V))-i{cosh(2/3t)((7V)- 



1) + [cosh2(2/ît)((iV) - 1)2 + 2{N) - (^)2]i/2}}. From 
this relation and Eq. (^2|) one can calculate the concur- 
rence as a function of température and the mean fermion 
number (N). This function is represented in Fig.2. We 
observe that the entanglement becomes maximal when 
the mean fermion number is 1 for fixed température T. 
Let us first comment on the n i— > 1 — n, or in view of 
Eq.(p2[) equivalently p i— > —p, symmetry of the function 
C. This fact can be understood from Eq.(|Ï2[). Clearly 
the only thing to check is that \z{t, p)\ — \z{t, —p)\. This 
latter statement easily follows by the use of particle- 



hole transformation i.e., c, <-!■ c 
z{t,-p) 



(j = 1, . . . , L) reahzing 



J 7 

z{—t,—p) along with the unitary mapping 



Cj (^l)"* Cj that, for bi-partite lattices, changes the 
sign of t. Another important feature is the existence of a 
threshold température, after which the entanglement dis- 
appears. Remarkably this threshold température is inde- 
pendent of the mean fermion number. This phenomenon 
is related to the independence on an external mag netic 
field displayed by the associated spin model |p^ , |ï^ . 

For large L, we plotted in Fig. 3 the concurrence as a 
function of température for différent p. We observe that 
the threshold température is independent of p. Moreover 
it is worth noticing that, for sufRciently high p ,i.e, filling, 
one has a non monotonie behavior of the concurrence as 
a function of T. Indeed we see that the entanglement 
can initially increase as the température is raised. This 
phenomenon is due to the fact that the chemical poten- 
tial in fermionic Systems plays a rôle analogous to the 
external magnetic field Bz for spin Systems When 
Bz is large enough the ground state is given by a product 
state (ail the spins aligned), hence entanglement in the 
thermal state is due to the excited eigenstates. Of course 
for T large enough one always gets C = in t hat the 
Gibbs state is approaching the maximally mixed state. 



IV. ENTANGLEMENT AND 
SUPERCONDUCTIVITY 

In this section we discuss fermionic entanglement in 
simple superconducting Systems and explore the relations 
between entanglement and superconducting corrélations. 
Let us first consider a BCS-like i.e., mean- field, model 



A. BCS-like superconductivity 

The foUowing BCS-like Hamiltonian describes the pair- 
ing between fermions carrying momentum k and —k, 

H ^ J2kHk, 

Hk = Êfc {rik + n^k - 1) + Afc c\c^_f, + CfcC_fc. (24) 

The quantifies = |Afc|e*'^'= are order parameters of 
conductor-superconductor phase transitions. They are 
determined by the self-consistent relations Aj, = (ckc^k)- 
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and above a critical température they vanish thus signal- 
ing the absence of superconducting corrélation. 

The structure of Eq. ( p^ ) clearly suggest that the rele- 
vant tensor-product structure of this problem is given by 
'Hf — ®k{hk ® /i-fc), where := span{(c|,)" |0) / a = 
0, 1}. Moreover it is very simple to check that, for any fc, 

the operators Jk := Uk + ri-k - 1, J'^ ■= c[-'^-fe' 
C-kCk span a su{2) Lie-algebra. The states \aa) :— 
\a)k <E) \a)-k (o = 0, 1), realize a spin-1/2 représentation 
of such an algebra. Therefore the Hamiltonian (|2^), that 
it is équivalent to a spin-1/2 particle in an external mag- 
netic field along the x direction, is readily diagonalized. 
For example, if we define 0^ := tan^^ |Afc|/efc, the ground 
State is given by (^k\—)k where 

cos(0fc/2) |00)fc - e'*'' sin(0fc/2) |ll)fe. (25) 

It corresponds to the eigenvalue Ek- = —Ek- Here E^ — 
^/e^+\Ak\^. The self-consistent équation for the A^'s 
reads 



sinh(/3gfc)sin(gfc) 
2[cosh{f3Ek) + 1] 



(26) 



The concurrence associated to the thermal state Pk{f3) 
exp{-(3 Hk)/Z is given by 



C 



max{0, Sïnh{P Ek) sm{6k) — 1} 
œsh(JËk)Tl 



(27) 



In the limit of T ^ 0, the concurrence becomes sin(6'fe) 
which is just the concurrence of the ground state |— )fc. 
By solving thèse équations the numerical results are given 
in Fig.4. One can see that the concurrence goes to zéro at 
a température slightly lower than superconductor critical 
one. Since in the température range [0, Te] the function 
Afc(T) is invertible one can express the concurrence as a 
function of the order parameter only. This is illustrated 
in Fig. 5 which shows that it is necessary to have a certain 
amount of superconducting corrélation in order to have 
an pairwise entangled thermal state. Notice that we set 
efe = that in a grand-canonical picturc {ek ^ tk ~ pi) 
corresponds to half-filling 

We would like to observe now that a non mean-field 
Hamiltonian formally analogous to (|4|) plays an im- 
portant rôle in the excitonic proposai for QIP by Bi- 
olatti et al pï] . In that case the fermionic bilinear 
terms cl^. are replaced by Xi := cldl^{i := {l,m)) 

where c| (dl^) créâtes an électron (hole) in the ^-th (m- 
th) state of the conduction (valence) band of a semi- 
conductor. The order parameter A^ becomes a (in- 
dependently controUable) coupling to an external laser 
field. The excitonic index l can be associated to L dif- 
férent, spatially separated, quantum dots; this implies 
that l ^ V ^ [Xi, Xii] = 0. Let |0) dénotes the particle- 
hole vacuum (ground state of the semiconductor crystal) . 
Since Xf = one immediately sees that the "excitonic 
Fock space" span{J|; X"' |0) / ni — 0, 1} is isomorphic 



to a L qubit space |2Ï|]. This example shows that one 
can consider spaces allowing for a varying number of 
"particles" that nevertheless are fuUy legitimate quan- 
tum state-spaces. No super-selection rule violation (pos- 
sibly due to spontaneous symmetry-breaking) has to be 
invoked. Next we consider another kind of , non mean- 
field superconductivity, i.e., the 77-pair superconductivity. 



B. 7;-pair superconductivity 

Yang 1^ discovered a class of eigenstates of the Hub- 
bard model which have the property of off-diagonal long- 
range order (ODLRO) which in turn implies the 

Meissner effect and flux quantization |^^. Let us begin 
by introducing the r/-operators 



3 = 1 



L 
2 



L 

i=i 



(28) 



which form the su(2) algebra and satisfies [rj, 77^] = 
2r]^ , [77^,77^] = ±rj^ . Here the fermions have spins and 
the operator nj — c'j^Cji -|-c^|Cj|. The operators 77^ also 

satisfy the relations (t;*)^^^ = 0, which reflect the Pauli 
principle, i.e., the impossibility of occupying a given site 
by more than one pair c^j^c||. 

In this context the relevant state-space is the one 
spanned by the 2^ basis vectors 



\ni, = l[{clclr\yâc){ni = 0,1 yi). (29) 

1=1 

From the above basis it is évident that the space 
span{|7ii, 77i)} is isomorphic to the L-qubit space. 
This is easily seen by defining the mapping 

A' := I7.1, ...,nL) ^ ®ti|7.i) = ®ti(a+)"'|0). (30) 

Then we can produce L 'number' state | A'^) by applying 
successive powers of 77+ on the vacuum state defined by 
77IO) = 0. So 



\N) ^Af-'/^f]+^\0), N^O,. 



(31) 



whereA/' = L!A^!/(L-iV)!. The span ofthe |iV)'s {N = 
0, . . . , L), known as 77-paired states, forms an irreducible 
spin-L/2 représentation space. A mentioned above what 
makes the number states interesting is the fact that t hey 
have been shown to have ODLRO. Indeed from Eq.(|3l|) 
the foUowing, distance independent, corrélation function 
is obtained 

O^.L = {N\c],c]^ci^cn\N) = (32) 

In the thermodynamical limit (N, i — > 00 with N/L — n 
) C^N.L goes like to n(l — n), which is nonvanishing as 
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|j — Z| oo. In other words the number state exhinits 
ODLRO, and thus is superconducting. 

The two-site reduced density matrix has the form ( |ÏÔ| ) 
in which 



{L - N){L - N - 1) 



O 



N,L, 



L{L-1) 

_ N{N- 1) 



L(L-l) ■ 

From thèse relations and Eq. (|2|) one finds 

{N-l){L-N-iy 



C = 2{ On,l - 



N,L- 



L{L-l) 



112- 



(33) 



(34) 



Notice that the formula above could have been directly 
obtained from Ref ( jî^]) where the entanglement be- 
tween any pair of qubits in a Dicke state has been com- 
puted. Indeed by means of the the mapping (pÔ|), one 
can identify the number state with the usual Dicke state 
In the thermodynamical limit one has u i— > (1 — 
n)^, V i— > T? , z = Oiv.L ^ niX — n) thus, from Eq. jï^ ) 
the entanglement becomes zéro. So we see that the pair- 
wise quantum entanglement does not exists although we 
have 77-pairing superconductivity in the number state. 
We finally notice that in the 77-pair cohérent states dis- 
cussed in |Q, there is ODLRO, but being a products 
there obviously no entanglement. 



V. CONCLUSIONS 

The Fock space of many local fermionic modes can be 
mapped isomorphically onto a many-qubit space. Using 
such a mapping we studied entanglement between pairs of 
(spinless) fermionic modes. This has been done both for 
the eigenstates and the thermal state for a model of free 
fermions hopping in a lattice (entanglement between local 
modes as function of température and filling). In the 
free fermionic model we analysed entanglement between 
local modes as function of température and filling. In 
particular we found that above a threshold température 
the thermal state becomes separable. 

We studied the relations between pairwise entangle- 
ment and the superconducting corrélation in both the 
BCS-like model and the 77-pair superconductivity. For 
the BCS-model, that a finite value of the superconducting 
order parameter is required to obtain entanglement in the 
thermal state. Notice that this last statement establishes 
a direct connection bewteen quantum entanglement and 
a real phase transition |^5| . For ?7-pair superconductivity 
we found that pairwise entanglement is not a necessary 
condition for the superconductivity. 

Despite their simplicity, our results seem to suggest 
that the quantum information-theoretic relevant notion 
of quantum entanglement can provide useful physical in- 
sights in the physics of many-body Systems of indistin- 
guishable particles. 
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TABLE L The concurrence for three spécial values of mean 
fermion number. 
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FIG. 2. The concurrence as a function of the température 
and the mean fermion number for L = 2 and t = 1. 
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FIG. 1. The concurrence as a function of the filling n in an 
infinité lattice 
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FIG. 3. The concurrence as a function of the température 
for différent jj.: jj. — 0.1 (cross points) , fi = 1.0 (circle points), 
and /i — 2.0 (diamond points). The parameter L = 100. 



7 




0.05 0.1 0.15 0.2 0.25 



FIG. 4. The order parameter (cross points) and the con- 
currence (box points) as a function of the température. The 
parameter ek = 0. 
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FIG. 5. The concurrence against the order parameter. The 
parameter €k = 0. 



